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Abstract 

We discuss issues related to the resummation of non-global observables in 
QCD, those that are sensitive to radiation in only a part of phase space. 
Examples of such observables are certain single-hemisphere event shapes in 
e + e~ and DIS. Compared to global observables (those sensitive to all emis- 
sions, e.g. the e + e~ thrust) a new class of single-logarithmic terms arises. 
These have been neglected in recent calculations in the literature. For a 
whole set of single hemisphere e + e~ and DIS event shapes, we analytically 
evaluate the first such term, at order a|, and give numerical results for the 
resummation of these terms in the large- Nq limit. 
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1 Introduction 

It has long been known that in certain restricted regions of phase space the coefficients 
of the QCD perturbation expansion can be parametrically large. For semi-inclusive 
observables this is often due to the incomplete cancellation of soft and/or collinear 
logarithms. 

A set of variables particularly sensitive to this kind of problem is event-shape distri- 
butions which we shall use as an illustration throughout. Typically one finds that 
the perturbative series is dominated by terms (a™ ln 2n_1 v)/v, where d C 1 denotes the 
value of the event-shape variable. In order to have a meaningful perturbative descrip- 
tion these terms have to be resummed to all orders. This has been done for a range of 
observables 

Most of the observables considered so far have the property of exponentiation, im- 
plying the following structure for the resummation, 



where the distribution is given by dH/dv. The remainder function D(v) has the property 
that it goes to zero for v — > 0; L = \nl/v and a s = a s /2ir. The terms resummed by 
Lgi(a s L) are referred to as leading (or double) logarithms (LL), while those in g2(a s L) 
are next-to-leading (or single) logarithms (NLL). For certain variables, rather than a 
single exponent one finds a sum or integral over exponents ]7} [| . 

A property common to most of the observables considered until recently is that they 
are sensitive to emissions anywhere in phase space. We call these global event shapes. 
Examples are the thrust, the heavy-jet mass, the C-parameter. In this paper we are 
interested in non-global observables, those sensitive to radiation in only a part of phase 
space. We shall discuss mainly single-hemisphere event-shapes in e + e~ (especially the 
jet-mass). Other examples include a number of DIS current-hemisphere observables 
and quantities such as the mass and broadening of single jets in multi-jet ensembles. 

The literature on the subject of global event shapes often makes use of a combination 
of results in separate hemispheres to obtain the properties of the whole event |], p], 
[5|. When the separate hemisphere results are combined in this fashion the answer 
obtained is correct. However this can be misleading, since it turns out that the existing 
single-hemisphere results, applied literally to just a single hemisphere event-shape, are 
incomplete at NLL level. This has not been taken into account in the recently presented 
results for the light-jet mass and the narrow-jet broadening |7j and for the narrow-jet 
thrust-minor M . In the case of the first two variables the incompleteness of the results 
can be seen from a comparison of the O (a 2 L 2 ) term of the expansion of the resummation 
with fixed-order predictions from a Monte Carlo program such as Event2 ]TU| 

2 Here we refer to the original version of [jTj where it was stated that terms down to a™L 2n ~ 2 of the 
expanded result were under control. In a revised version (v2) which appeared on the hep-ph arXiv 
shortly before the submission of our paper, the claims regarding the accuracy were revised to state that 
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To illustrate the point it is convenient to study a simpler variable, the hemisphere 
jet-mass, defined as 



(2) 

where the hemisphere is taken with respect to the thrust axis. Once the distribution of p 
has been determined, it can be combined with the known heavy-jet mass distribution Q 
to obtain the light-jet mass distribution, using the following exact relation 

da da da da da 

+ 1 = 1 + 1 = 2 ^' ( 3 ) 



^Plight ^Pheavy ^Pleft fright dp ' 

where p is the hemisphere jet-mass. An analogous relation applies for the jet-broadenings. 
For the narrow-jet thrust minor the situation is more complex but considerations similar 
to those addressed here will apply. 

One can understand the origin of the qualitative difference between the heavy-jet 
and single-hemisphere mass resummations by examining figure |], which depicts two 
soft large-angle gluons, whose energies are ordered Q 3> fcj 3> k 2 - For the heavy-jet 
mass there is a cancelation between the two contributions because the emission of k 2 
does not affect the value of pheavy (which is determined by the harder emission, k\). 
In the case of the right-jet mass, k\ has no effect on p r i g ht, so the event shape receives 
a contribution only from the softer emission, k 2 - This spoils the cancelation with the 
virtual diagram, figure [L]a, thus leading to a term a 2 L 2 , i.e. a NLL effect. Such effects 
are a general feature of non-global observables. 

At all orders, to NLL accuracy, the form of the resummed jet-mass distribution will 

be 

£(p) = (l + a s C[ q) ^ S(a s L) S 9 (a s , L) + a s C[ g) E g (a s , L) + D(p) . (4) 
£„ and E„ are the resummed jet-mass cross sections which can be obtained by integrat- 



ing the quantities J q and J g first derived in |L1] and widely used in the literature. They 



have structures like that of the exponential factor in eq. ([!]) . The sum of two resummed 
contributions arises because in a fraction O (a s ) of events the hard particle in the right 
hemisphere (7Yr) is a gluon. 

The new effect considered in this article is embodied by the function <S.Q We write 
its expansion as 

S(a s L) = l + J2Sn(a s L) n . (5) 



only terms down to a™L 2 ™ -1 are controlled. In our notation that corresponds to controlling only the 
a s L term of g2, but not terms a"L n with n > 2, nor C\. Here we follow the convention adopted in pj] 
rather than that in and use the term next-to- leading to refer to logarithms in the exponent (s), not 
in the expansion of the exponent. 

3 In principle a similar function multiplies E fl , but at NLL accuracy it can be neglected. 
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It is straightforward to exactly compute the first non-trivial term S2 and this is done 
in the following section. The full computation of S involves considering an ensemble of 
an arbitrary number of large-angle energy-ordered soft gluons in fix,, which coherently 
emit a single, softer gluon into 7Yr. For reasons elucidated later it is difficult to carry 
out an all-orders treatment of such an effect analytically. We therefore opt to treat these 
effects using a Monte Carlo algorithm valid in the large- Nc limit. This is outlined in 
section ^| and further details are given in the appendix. 

Finally in section El we compare our results to the O (a?) predictions of Event2. 



Phenomenological predictions including this effect will be shown elsewhere [12 



2 Fixed order calculation 

First we calculate the contribution to the jet-mass distribution from the configuration 
in figure ^Jo, considering the right-hemisphere jet for concreteness. We introduce the 
following particle four-momenta 

fc a = §(l, 0,0,1), (6a) 

fc 6 = |(l, 0,0,-1), (6b) 

ki = £i^-(l, 0, sin #1, cos 6\) , (6c) 

^2 = #2"^(lj sin ^2 sin 0, sin 9 2 cos0, cos 62) , (6d) 

where we have labelled the quark and antiquark as a and b and defined energy fractions 
%i,2 <C 1 for the two gluons. We have ignored recoil in the kinematics, because the 
jet-mass is insensitive to it. 

When gluon 2 is in the jet mass has the value p = 22(1 — cos#2)/2. When only 
the quark is in 7Yr, p = 0. 

We write the matrix element for ordered two-gluon emission as (see for example [ |i~3|j ) 
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(al)(16) V 2 (a2)(21) 2 (62) (21) V 2 / ( a2 )(26)y v ; 

= C 2 + C F C A W 2 , (7b) 

where (ij) = h L ■ kj. The result is valid for 1 3> £i 3> £ 2 as well as for the opposite 
ordering of the gluons, and furthermore is completely symmetric under interchange of 
k\ and k 2 . We have however chosen to write it in an asymmetric form to emphasise 
the dipole structure of the emissions, namely radiation of gluon k\ from the ab dipole, 
followed by the radiation of gluon k 2 from the al, lb and ab dipoles. 

The Cp piece of the matrix element, W\ corresponds to independent gluon emission 
and is included in the usual resummation of the quark jet-mass, i.e. in E g of eq. (H). To 
study specifically the modification relative to the standard quark-jet mass result due to 
configurations like that shown in figure |I]b, one must consider the CfCa part, W 2 , of 
the emission probability above. 

The integral to be considered, related to the probability for the jet mass to be less 
than p, is then 



/a s \ 2 f° f r ' 2 " ■ r ' 2r ' 

— CfCa J ^ c °s^i J dcos9 2 



2ir L 2ir 



Q 4 f 1 f 1 fx 2 \ 

— J x 2 dx 2 J x x dx x 6 - cos 9 2 ) - pj W 2 , (8) 

whose dominant term at small p is S2 «g In 2 1/p. 

The integrals over the energy fractions are straightforward. Keeping only the In 2 1/p 
piece and performing the azimuthal average we have 



,1 ,0 

S 2 = —2CfCa \ dcos9 2 / dcos^! fi 2 (cos^i, cos# 2 ) 



(9) 



where the angular function fl 2 is 
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fi 2 = • (10) 

(cos# 2 — cos#i)(l — cos^i)(l + cos# 2 ) 



Integrating over the polar angles we obtain 



TV 2 



S 2 = -C F C A -. (11) 



We note that possibly related u ir terms have been discussed in [14 



3 Resummation 

To perform the all orders resummation of S(a s L) we need to consider general configu- 
rations of soft gluons. We define the probability Pc{L) for a given configuration C at a 
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resolution scale L {i.e. not resolving gluons with energies less than Qe~ L ). Given Pc, 
the probability Pc of a configuration C with one extra gluon with scale L' > L, polar 
angle 9' and azimuth <$' is (see for example |L5] ) 

Pc>(L') =a s {L')A c (L,L , )Fc(6 , ,<t>')Pc(L), (12) 

where the form factor is given by 

lnA c (L,L') = -J dL" J dcosO d(j) a s (L") F e (6,(j)) , (13) 

and Fc{6,4>) describes the angular and colour-structure of the radiation pattern from 
configuration C. The angular integrations in (|13D diverge whenever the angle is collinear 
to one of the emitting particles. In practice it is therefore convenient to introduce an 
angular cutoff e to regulate these divergences both in the real emissions and the virtual 
corrections. 

We obtain S by calculating the probability of there being no emissions in 7Yr down 
to a scale L, divided by the corresponding probability had the only source of emissions 
been the original qq pair. This gives 

S(a s L) = ^=L== £ ML), (14) 

V A «H L i C|«r empty 

where the sum runs over all configurations C which contain no emissions in 7Yr. 



3.1 Large- AT C limit and Monte Carlo implementation 

In practice two considerations make it difficult to implement the above approach ana- 
lytically. One is that the colour algebra involved in the determination of Fc becomes 
progressively more complicated as the number of gluons increases. The other problem 
is simply that the treatment of the geometry of the many large-angle gluon ensem- 
ble quickly becomes prohibitive. The first problem can be partially solved by taking 
the large- Nc limit. To address the issue of the geometry we shall use a Monte Carlo 
approach. 

In the large- Nc limit, one can represent gluons by pairs of colour- anticolour lines, 
as illustrated in figure When squaring the amplitude one ignores contributions that 
in terms of their colour flow are topologically non-planar, because they are suppressed 
by powers of 1/N C |Tj|. As a result, for an ensemble of n gluons, Fc just reduces to a 



sum of independent emission intensities from n + 1 independent dipoles: 

2C A 



E 



.. (1 - cos0 ife )(l - COS0 kj ) 

dipoles— i] J 



(15) 



When the dipole ij radiates a gluon k it splits into two dipoles, ik and kj. Thus the 
dipole structure is determined by the history of the gluon branching. 
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Figure 2: Left: the kind of diagram which must be considered in the calculation of S. 
Right: the same diagram represented in the large- Nc limit, with gluons shown as pairs 
of colour lines and quarks as single colour lines. 



Such a branching pattern can be very naturally implemented using a Monte Carlo 
algorithm. At first sight one might envisage calculating the two factors in sepa- 
rately, using the Fc to generate the distribution of radiation for each new configuration. 
However because of the collinear divergence along the direction of the quark in T-Cr, 
only a tiny fraction of events would be free of emissions in TIr and so contribute to the 



sum in ([TJ|) . The sum would therefore have a large relative error, which would translate 
to a large absolute error on S because of the division by the small quantity a/ A a b(L). 

Instead a more efficient procedure involves moving the division by a/ A ab (L) directly 
into the calculation of the Pq- This can be achieved using a modified radiation intensity, 
Fc (for both the emissions and the virtual corrections), 



F ab (9,<P)e(9) 



(16) 



where one subtracts out the radiation intensity F a b which would have been produced by 
the original qq pair (in the large- Nc limit). One calculates quantities Pq using analogs 
of eqs. fll2|) and flT3|) with Fc replaced by Fc and then S is simply given by 



S{a s L) 



C\Hn empty 



(17) 



It should be kept in mind that since Fc is negative in certain regions of phase space 
one loses a strict probabilistic interpretation for the Pc- Nevertheless the sum over 
configurations is well-defined and meaningful. 

The exact details of the Monte Carlo algorithm are given in the appendix. Here we 
restrict ourselves to giving a parameterisation for S obtained by fitting to the Monte 
Carlo results: 



S(a s L) ~ exp -C F C A — 



7T" 



(at)' 



with 



1 f 1 dx 

t(a s L) = — — a s (xQ) 

Z7T J e -L X 



1 



(bty 



In 



r 



i 



4tt/3 1 - 2f3 a s L 



(19) 
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where (3 = (11C A - 2n / )/(12vr) and 



a = OMC A , b = 0.86C A , c= 1.33. (20) 

The parameterisation should be accurate to the order of a few percent (better in most 
of the region) for t < 0.7, corresponding to 1 — 2a s p L > 0.005.0 

Actually, for the purposes of the fit one replaces CpC A in ([18]) with C\/2 since the 
Monte Carlo works in the large- Nc limit. But for use in phenomenology one wishes to 
have the exact colour structure at least at O (a 2 ), hence the use of CpC A in (|T8|). 



4 Checks and conclusions 

It is useful to check our results against fixed order results from the next-to-leading order 
Monte Carlo program Event2 [ 10(1 . First it is necessary to determine the constant terms 



C[ gS> and C[ 9 \ which are obtained by requiring consistency between ([|) and a full O (a s ) 
calculation. It is straightforward to show that they are given by 

C[ q) = \ {CI - r 3 ) , C[ 9) = | , (21) 

where C\ is the constant determined for the thrust in JI]] and a s r 3 is the probability of 
there being a hard gluon along the thrust axis, as given in :f] 
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it' 2 


4) 




+ T 





(22) 



The plots in figure |3| show the difference between the O (a 2 ) exact and resummed 
distributions. To NLL accuracy we aim to account for terms down to a 2 L 2 in S and 
therefore a 2 L in the distribution. Accordingly the difference should be independent of 
p for small p. One can see that this is the case for each individual colour factor. We 
also show that without the contribution S calculated in this paper the CpC A part of 
the result is not consistent with Event2. 

Using eq. (|3|) it is possible to combine our results for the single-hemisphere mass 
together with those of |lj for the heavy-jet mass in order to obtain predictions correct 
to single- log level for the light-jet mass (including the appropriate C\ terms). 

We point out that the form computed for S applies to all single-hemisphere 2-jet 
event-shapes that for large-angle gluons have a value of the order of the transverse 
momentum: in e + e~ this means the hemisphere broadening (for the constant terms a 
relation analogous to (^T|) applies, with C\ replaced with C*f T ). In DIS our form for 
S applies to the current hemisphere jet-mass, the C-parameter, and the thrust and 

4 The accessible range of t is limited by two issues: firstly only a small fraction of events are generated 
at large t, requiring considerable statistics in order to investigate that region; and secondly because an 
accurate determination of S at large t requires a very small angular cutoff, which leads to there being 
many dipoles in an event, and a consequent slowing down of the evolution. 

D Only in the hep-ph arXiv version 1. 
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Figure 3: The difference between fixed order results from Event2 and the expansion of 
the resummed results. The left-hand plot shows results for the coefficient of o£,CfCa 
- the two sets of points correspond to neglecting or accounting for S. The right-hand 
plot shows the coefficients of the CfTj and Cp pieces. 



broadening with respect to the thrust axis. The detailed phenomenology of the DIS 



variables will be considered elsewhere ||12|| . It indicates that the inclusion of S reduces 
the peak height (in da /dp) by about 30% before matching to the O (af) calculation, 
and by about 20% after matching (using an i?-type matching [l|]). 

For the narrow-jet K out and jet masses in multi-jet events (including hadron-hadron 
event shapes) different forms will apply because of the more complicated geometry 
of the underlying hard event, but the Monte Carlo algorithm described here can be 
adapted to those well. 
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Appendix 



Below we give details of the Monte Carlo algorithm used to calculate S.f\ 

For each dipole T>i, the collinear angular cutoff e introduced in section |3] defines a 
rapidity interval Arji for emission in the dipole c.o.m. frame. We shall denote by T>\ 
the dipole which has as one of its sources the quark in Ti^. The evolution variable used 
will be t, defined in (19|). This will account for the running of the coupling. 

One starts an event with a quark- ant iquark dipole (labelled ab), t = and a weight 
w = 1. One then goes through the following steps to determine an event's contribution 
to the quantity —dS/dt: 



1. Establish A?7 tot = Yli^Vi- Increase t by an amount At chosen according to a 
distribution proportional to 

exp {-2C A A Vtot At) . (23) 



2. Choose a dipole at random, such that the probability of obtaining T>i is Arji/Arjtot- 
Create a gluon uniformally in rapidity and azimuth in the dipole centre of mass 
frame, such that in the lab frame its angle relative to either of the dipole sources 
is larger than e. 

3. There are now various possibilities: 

— The gluon is in Hl- insert it into dipole i so as to create two new dipoles 
(removing the old one from the list). Go to step 1. 

— The gluon is in Ti^ and was not emitted from T>\. add w to the bin corre- 
sponding to the value of t and start a new event. 

— The gluon is in 7Yr and was emitted from T)\. denote the momenta of its 
sources by pi and p a (the latter is just the 7Yr quark momentum) and the 
new gluon momentum by p g . One determines the value of 

x = i _^ = i 

Fat (gb) (a£) • 

If X > then with probability X we add w to current t-bin and start a new 
event; with probability 1 — X we throw away the gluon that has just been 
generated and go to step 1. If X < then we add Xw to the current t bin, 
multiply w by 1 — X, and go to step 1. This procedure is a simple extension 
to negative probabilities of the veto algorithm and corresponds to replacing 
Fc, used implicitly in the previous steps, with Fq. 

6 We present a version of the algorithm which has weighted events; there exists also a slightly more 
efficient version of the algorithm with unweighted events, but it requires some integrations to be done 
analytically and is therefore less easily generalisable to observables with more complex boundaries in 
phase space such as jet masses in multi-jet events. 
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